#include<bits/stdc++.h>
using namespace std;
namespace IO{
    char buff[1<<21],*p1=buff,*p2=buff;
    char getch(){
        return p1==p2&&(p2=((p1=buff)+fread(buff,1,1<<21,stdin)),p1==p2)?EOF:*p1++;
    }
    template<typename T>
    void read(T &x){
        char ch=getch();int fl=1;x=0;
        while(ch>'9'||ch<'0'){if(ch=='-')fl=-1;ch=getch();}
        while(ch<='9'&&ch>='0'){x=x*10+ch-48;ch=getch();}
        x*=fl;
    }
    template<typename T,typename ...Args>
    void read(T &x,Args& ...args){
        read(x);read(args...);
    }
    char obuf[1<<21],*p3=obuf;
    void putch(char ch){
        if(p3-obuf<(1<<21))*p3++=ch;
        else fwrite(obuf,p3-obuf,1,stdout),p3=obuf,*p3++=ch;
    }
    char ch[100];
    template<typename T>
    void write(T x){
        if(!x)return putch('0');
        if(x<0)putch('-'),x*=-1;
        int top=0;
        while(x)ch[++top]=x%10+48,x/=10;
        while(top)putch(ch[top]),top--;
    }
    template<typename T,typename ...Args>
    void write(T x,Args ...args){
        write(x);write(args...);
    }
    void flush(){fwrite(obuf,p3-obuf,1,stdout);}
}
using namespace IO;
const int N=1e3+5;
bool vis[N];
int p[N],tot;
void init(int x){
	for(int i=2;i<=x;i++){
		if(!vis[i])p[++tot]=i;
		for(int j=i+i;j<=x;j+=i)vis[j]=1;
	}
}
int n,mod;
inline int Mod(int x){return x>=mod?x-mod:x;}
int a[N][N],b[N][N],c[N][N];
int A[N],B[N],C[N];
void Init(int x){
	int sum=1;
	for(int i=1;i<=n;i++){
		sum=1ll*sum*p[x]%mod;
		A[i]=B[i]=sum;	
	}
}
void Work(){
	fill(C+1,C+n+1,0);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			C[j]=Mod(C[j]+1ll*A[i]*a[i][j]%mod);
	for(int i=1;i<=n;i++)A[i]=C[i];
	fill(C+1,C+n+1,0);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			C[j]=Mod(C[j]+1ll*A[i]*b[i][j]%mod);
	for(int i=1;i<=n;i++)A[i]=C[i];
	fill(C+1,C+n+1,0);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			C[j]=Mod(C[j]+1ll*B[i]*c[i][j]%mod);
	for(int i=1;i<=n;i++)B[i]=C[i];
}
bool check(int x){
	Init(x);
	Work();
	for(int i=1;i<=n;i++)
		if(A[i]!=B[i])return 0;	
	return 1;
}
inline void sol(){
	read(mod,n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			read(a[i][j]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			read(b[i][j]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			read(c[i][j]);
	int flag=1;
	for(int i=1;i<=tot;i++){
		if(!check(i)){
			flag=0;
			break;
		}
	}
	if(flag)putch('Y'),putch('e'),putch('s');
	else putch('N'),putch('o');
	putch('\n');
	return;
}
int main(){
	freopen("matrix.in","r",stdin);
	freopen("matrix.out","w",stdout); 
	init(100);
	int T;
	read(T);
	while(T--)sol();
	flush();
	return 0;
}

